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In a diffusive conductor the eigenstates are spread over the entire sample. However, with certain 
probability, an anomalously localized state (ALS) can occur, i.e. the wave function assumes anoma- 
lously large values in some region of space. Existing analytical theories of ALS are based on models 
described by a continuous (Gaussian) random potential. In the present paper we study ALS in a 
■ lattice (Anderson) model. We demonstrate that close to the center of the band, E = 0, a new type 

of ALS exist and calculate analytically their likelihood. These ALS are lattice-specific and have no 
analog in the continuum. Our findings are relevant to numerical simulations, which are necessarily 
performed on a lattice. We demonstrate that inconsistencies with "continuous" results reported in 
the previous numerical work on ALS can be explained within our analytical theory. Finally, we point 
out that, in order to compare the numerics with the "continuous" ALS theories, simulations must 
be carried out not too far from the band edges, within the band, where the continuous description 
applies. Simulations performed for E close to the band center reveal lattice-specific ALS that do 
not exist in continuous models. 
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Introduction. In a weakly disordered conductor the typical value of an eigenfunction intensity, (r) | 2 , is of order 
L~ d , where L is the sample size, in d dimensions (d = 2, 3). However, with certain probability, the intensity can assume 
anomalously large values. The study of such rare events in diffusive conductors was pioneered in Ref. [1] and further 
pursued in Refs. [2-4]. The " prelocalized" states, studied in Refs. [1-4], exhibit an anomalous buildup of intensity in 
some region of space, of a size larger than the mean free path, I. The properties of these states are universal, in the 
sense that the disorder enters only via the mean free path. Another type of rare events was identified and studied 
in [5]. The corresponding eigenstates, designated as "almost localized states", are confined primarily to small rings, 
of a sub-mean-free-path size. These states are non-universal, i.e. sensitive to the microscopic details of the system. 
\q . In particular, their likelihood sharply increases with the correlation radius, R c , of the disordered potential (for fixed 
' value of I). In what follows we designate any type of an anomalously large buildup of intensity as an anomalously 
, localized state. 

The above mentioned analytical studies of the ALS were limited to models described by a continuous (Gaussian) 
random potential. On the other hand, numerical studies of disordered electronic systems are necessarily performed 
' on the lattice, most often within the Anderson model [6] , with the tight-binding Hamiltonian 



H = 5^4tV+5^V r 4cp, (1) 



where cj. is the creation operator of an electron at site r of a ci-dimensional hypercubic lattice with lattice constant 
equal to 1, and V r is a random on-site energy with r.m.s. = (V 2 ) 1 / 2 . The Anderson model has become a 
powerful tool for numerical study of various disorder-related phenomena. As the computing capabilities constantly 
grow, allowing diagonalization of large-size matrices, more and more accurate information can be inferred from the 
simulations. As a result, the early success [7] in confirmation of the scaling theory [8] was followed by recent numerical 
studies that have successfully addressed more delicate issues such as (a) critical exponents [9] and critical behavior 
of the eigenfunctions in 3d [10], (b) quantitative characteristics of the quantum Hall transition [11], (c) different 



aspects of the level statistics at the Anderson transition [12,13], (d) Anderson transition in 2d [14] possible with 
spin-orbit coupling [15], (e) the critical conductance distribution at the transition [16,17], (f) verifying scaling for the 
full conductance distribution [18]. 

This successes have encouraged a number of authors [19-25] to employ the Anderson Hamiltonian for numerical 
study of the ALS in disordered conductors. In particular, the subject of interest is the function fd(E,t) defined as 

/,,(£./) = —( > /('- <l>.,lO)\ 2 L d )6(E-E a )\, (2) 




where \l/ a and E a are the eigenfunctions and eigenenergies of the Hamiltonian (1), respectively, and v{E) is the 
density of states. ALS are responsible for the large-i tail of the disorder-averaged distribution (2) of the eigenfunction 
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intensity at a given energy, E. They correspond to the anomalous buildup of certain eigenfunctions inside the volume 
L d [1-4]. Simulations performed have revealed a number of unexpected peculiarities in the likelihood of the almost 
localized states: (i) in 2d, the behavior f 2 (E,t) oc exp (— Cb In 2 t) which is in accord with theoretical prediction of 
Refs. [1-4] was obtained [19]. However, upon changing the disorder magnitude, A2, the constant Ce did not scale 
with the conductance g ~ A^~ 2 ; (ii) the magnetic field dependence of Ce in Ref. [19] turned out to be very weak. 
This is in contrast to the simulations [26] of the eigenfuctions intensity distribution, in which the model of the kicked 
rotator rather than the Anderson model was studied. While the simulations [26] have also revealed exp (-Cln 2 t) 
behavior, the coefficient C in the presence of the time-reversal symmetry was almost two times smaller than in the 
limit when this symmetry was completely broken; (iii) Simulations of Ref. [25] suggest that the likelihood of the 
almost localized states is non- universal. Namely, fd(E,t) depends on the correlation radius of the random on-site 
energies, V r , for a given conductance, g, which is, again, in conflict with theory; (iv) Simulations in 3d [22,23] indicate 
that the wave function intensities, \^(m x ,m y ,m z )\ 2 , of ALS in a cube with a side L, plotted in the lexicographic 
order m — > L 2 (m x — 1) + L(m y — 1) + m z , are structured; a typical wave function represents a system of well defined 
and almost even-spaced spikes, each spike extending over a certain narrow interval of m. In contrast, diffusive wave 
functions |\&(m)| 2 , plotted in the same way, do not exhibit any structure [22,23]. 

In the present paper we demonstrate that the above peculiarities stem from the fundamental difference between 
the ALS in continuum and on the lattice. This difference is most pronounced for energies close to the band center, 
E = 0, where all the simulations [19-25] were performed. 

One should realize that the point E = is rather special and that "continuous" theories can break down in the 
vicinity of this point, even if the mean free path is much larger than the lattice spacing. For instance, in the ld case, 
it has been recently explained in Ref. [27] why the single parameter scaling, well established for weakly disordered 
continuous models, is violated near E = 0. Consider the Schrodingcr equation (E(k) — E) Af. + ]P fc , Uk-k'^k' = 
for the Fourier components, A^, of wave function, where E(k) = 2 (cos 

kx cos /by) is the bare spectrum in 2d; the 
Fourier components, U p , of the on-site disorder, V r , are described by the following correlator (U p U q ) = A 2 J2g $p,q+G, 
where G are the vectors of the reciprocal lattice. The terms with G / are responsible for the difference between 
the lattice and continuum. Indeed, as it was demonstrated in Ref. [27], with regard to these terms, the bare spectrum 
E(k) = 2 cos A: in ld must be considered as a two-band spectrum E = ±2cos/c within the reduced Brillouin zone 
\k\ < 7r/2. In other words, the umklapp processes, resulting from the terms G / 0, are small in parameter 1/Z, 
which does not have an analog in continuum. Away from E = this parameter is diminished by a factor ( j _E7| Z) for 
\E\ 3> Similar "non-universal" phenomena near E = should persist also in lattice models in higher dimension. 
In 2d, the center of the band corresponds to the four lines k x ± k y = ±7r, which constitute a square (see Fig. 1). 
Due to umklapp processes at (0, ±7r) and (±7r, 0), this square is the "right" Brillouin zone, which accommodates the 
"upper" (E > 0) and the "lower" (E < 0) bands. The umklapp related processes arc now small in parameter l/l 2 , 
which is, again, specific to the lattice. 

Below we identify a new type of ALS, specific for the lattice, and calculate analytically their density. The analytical 
theory allows to account for all the peculiarities (i)-(iv) uncovered in the numerics. The main idea is that the intrinsic 
periodicity of the Anderson model offers a possibility to localize the electronic states by the periodic fluctuations, V r , 
with a period 2, which cause a "dimerization" of the underlying lattice. In ld, such Peierls-like fluctuations [28] create 
a gap in the spectrum of the tight-binding Hamiltonian (1). In order to "pin" the center of localization of the in-gap 
state to a certain lattice site, a n phase shift in the periodicity should occur at this site, by analogy to the topological 
solitons [29] . Our prime observation is that similar fluctuations (period-doubling along each axis accompanied by a n 
phase shift) are capable to localize tight-binding electron in higher dimensions, without formation of a gap. It is still 
convenient pedagogically to start from the ld case, and then generalize the theory to higher dimensions. 

ld case. In the absence of disorder, the bare density of states is v(E) = n^ 1 (4 — E 2 ) . As was explained above, 
the period-doubling fluctuation V(°) (n =/= 0) = Vosign (n) exp (inn) with n phase shift at n = creates a localized state 
with energy, E G [— Vo, Vo]. Indeed, upon introducing a vector \ n = |*(2n),*(2n-l)| the Schrodinger equation 
if if = E^ with on-site energies V(n) takes the form 

(E-V o- 3 -a 1 )x n = o--X n+1 +o- + X n - 1 , (3) 

where Oi are the Pauli matrices. The corresponding eigenvector X n * ex P (* 7rn ~ ^71 \n\) decays both to the left and 

to the right from the site n = [29] with a decrement 71 (Vo) = arcsinh^/ j (Vq — E 2 ). The actual position of the 

localized state within the gap [— Vo, Vo] is governed by the on-site energy at the origin, V(n = 0). Namely, V(n = 0) 
and E are related as V(n = 0) = E [1 + e~ 71 / sinh7i]. 
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For a Id interval of a length, L, the fluctuation V^°\n) with V^(0) = V(n = 0) results in the buildup, r = 
\^(0)\ 2 /\^(£/2)\ 2 = exp [27i(Vo)£], of the eigenfunction at n = 0. However, the statistical weight of the fluctuation 
V^(n) is zero. Random deviations of the on-site energies, V(n), from ±Vo give rise to a certain distribution, 
Pi,c(E,t), °f buildups, r. To find this distribution, we note that deviations of V^°\n) from ±Vo result in the 
fluctuations of the local decrement, 71. Then the expression for V\,c can be written as 



/>oo />0 

Vi X (E,t)= dV! dV 2 ...P(V 1 )...P(V c/2 ) 
Jo J -00 



S (t - e s ) = (6{t- e s )) Vn = U> h c (E,lnr) , (4) 



where S = 2^ n 7i(^n) = 2^ n arcsinhy | (V 2 — E 2 ), and P{Vi) is the distribution function of energy of the site i. 
To proceed further, we notice that the function V\,c satisfies the following recurrent relation 

Pi,£(lnr)= J dV 1 P(V 1 )P llC -i(hiT-2 11 (y 1 ) ). (5) 
It is easy to express the solution of Eq. (5) in terms of auxiliary function 

poo 

h(K,E)= dVP0/)exp[-2m7i(V)] . (6) 

J E 

Then we have 

A,£(lnr)=/ d K e* Klnr [Ii(k)] £/2 = e~^ (E ^ '. (7) 



The analytical expression for log-probability of the buildup, Ti y c(E,r), can be obtained in the domain of C where 
the the main contribution to the integral (6) comes from small n. Then we can use the expansion Ii(n) — 7i(0) — 
i (71) k — (7J) k 2 , where (71 (E)) and {^\{E)) are the average decrement 71 and 7^, respectively. Substituting this 
expansion into (7), we readily obtain 



2/i (0) (7i 2 )-(7i)' 
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where lnTi = \ (71) £/h(Q). The domain of applicability of Eq. (8) is A X C > ln(r/Ti) > Ai/I 1 / 2 . The origin of T x 
and of the second term in Eq. (8) is the "prefactor" in the functional integral (4). 



2d case. Our main finding is that, in 2d, with the bare density of states v[E) = n 2 + K ^ + J^j ^ , where 

K is the elliptic function of the first kind, so that v{E) cx In | | ~ 1 at E— >0, a straightforward extension of the Id 
approach applies. Namely, the period-doubling fluctuation V(n x , n y ) = V$[ sign(n x ) exp {iim x ) + sign(n y ) exp (iim y )\ 
creates an ALS without opening a gap in v{E). Indeed, with V(n x , n y ) — V(n x ) + V(n y ) being additive, the solution of 
the 2d tight-binding equation is multiplicative, i.e. ^{n x , n y ) oc exp [(iir/2 — 72) (n x + n y )\. The decrement 72(^0) = 



arcsinhy ^ (4V 2 — E 2 ) differs from 71 due to the fact that, in 2d, the energy E is the sum of energies of motion 

along x and y. Now, analogously to the Id case, we consider a square with a side, Lj\f2, as shown in Fig. 1, and 
introduce a distribution, V2,c{E,t), of the probability that the buildup from the perimeter to the center along each 
path exceeds r. The reasoning leading to recurrent relation between p2,c{E, lnr) = t c V2,c{E,t) and P 2 ,c~i goes 
as follows. The perimeter site, i, of the square, C, is connected to perimeter sites of the square, C — 1, with one 

(C~ 1) (C — 1) 

horizontal and one vertical link, as illustrated in Fig. 1. Denote with 17 and t^_ 1 the values of buildup from 
these two perimeter sites to the center. The evolution of r along the horizontal and vertical links can be expressed as 
ln-r^P = lnr} £ ^ + 2^^ , and lnr^ = lnr^ ^ + 272"'^, respectively. This leads to the following relation 

P2.C ({In TV, In 7-^}) = J ...J {dViP(Vi)}x 

P2X-1 ({lnr, - 2 7 ^' i) ,lnr i _ 1 - 2 7 <"' i >}) . (9) 
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To make Eq. (9) closed, we recall that the actual buildup, T{ is the minimal of the horizontal and vertical values, i.e. 
Tj = nim-JY/j^, T Vt i\. Taking this fact into account, the solution for V2,c has the form similar to Eq. (7) for the Id case 



P 2)£ (hlT) = 



in In r 



Mk)\ 



£/2 



C/2 



(10) 



where the function 12(a) is defined by Eq. (6) with 71 (V) replaced by 72 (V). Using the small-K expansion h(a) = 
^2(0) — i (72) k — (7I) k 2 , we arrive at the following 2d generalization of Eq. (8) 



ln 2 (r/T 2 )/ 2 2 (0) 



2/2(0) (7 2 2 >-(72> 5 



^|ln{J 2 (0,£)}|, 



(11) 



where lnT 2 = ± (72) C ~ A 2 £/7 2 (0). The remaining task is to express the intensity distribution (2) through 
the distribution of buildups, r. We consider only the 2d case. For a given sample size, L, and the fluctua- 
tion size, C, the values t and r are related via the normalization condition for $(n), which can be presented as 
tL~ 2 [(27 2 )~ 2 + t- 1 (L 2 - C 2 )} = 1. Taking into account that 27 2 £ = lnr, we obtain 



t \L 



C 



ln 2 r 



Expressing r from Eq. (12) and substituting into Eq. (11), we get 



InV 



(12) 



^2,C = 



C E 



where 



hit -In 1 



C F . 



(CsftX 
I Lint / 



(72) £ 



+ T \\n{I 2 (0,E)}\, 



2 2/2(0) ( 72 2 >- (72)^ 



(13) 



(14) 



Further steps depend on the relation between L and t. For L 3> \ftj (72) ~ \/tg the second logarithm in Eq. (13) 
containing can be neglected. Then we get f 2 (E, t) — 4min£jF2,£, where the factor 4 accounts for the four 

quadrants. Performing minimization, we obtain fe(E,t) — t/Ce^o. t, where 77 < 1 is a numerical factor. For the 



Gaussian distribution P(V) this factor is equal to 77 = 



1 



0.76. Optimal value of C satisfies the 



4(tt-2) In 2 

condition C <C Lhit/y/i, so that the above assumption is justified. This assumption implies that the normalization 
of ^(n) is determined by the region L — C, i.e. outside the fluctuation. In the opposite case, L <C y/tg, both terms in 
Eq. (12) have the same order. In this case, the smallest possible fluctuation size C ~ L\nt/y/t should be substituted 
into Eq. (13). This yields 



h = 



C E 1 



(72)/ 
2Vt 



In 2 t. 



(15) 



Since Ce ~ g and (72) 
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-1/2 



it is easy to see that both L-dependent terms in Eq. (15) are small compared to 1. 



Thus, we again arrive at / 2 = Ce In 2 t with Ce given by Eq. (14). 

We now return to the peculiarities in the numerical results listed in the Introduction, and discuss them in light of 
the picture of ALS based on the period-doubling fluctuations, (i) The dependence of Ce on the disorder strength, 
calculated from Eq. (11) for Gaussian distribution with a standard notation for the r.m.s., A 2 = (V 2 ) 1 ^ 2 = Y2r x l 2 W , 
is shown in Fig. 2(a). For the range of 2 < W < 4, studied in Ref. [19], the slope of lnC_E versus lnIF~ 2 ~ In g varies 
within the range 0.67-0.83 for the energy interval \E\ < 2. As seen in the inset, the scaling Ce oc g is recovered from 
(11) at W < 0.5. (ii) Insensitivity of Ce to the magnetic field: for the period-doubling fluctuation along both axes, 
V(n x ,n y ), considered above, this insensitivity is an immediate consequence of the fact, that the magnetic phases in 
hopping matrix elements can be formally "absorbed" into the phase factors in the on-site values of the eigenfunction 
^(n x ,n y ). The underlying reason for such gauging out is that the fluctuation V(n x ,n y ) is separable, (iii) Sensitivity 
of fd(E, t) to the correlation of the disorder for a given conductance: it is obvious qualitatively that for the correlation 
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radius, R c , exceeding the lattice constant, the likelihood of the period-doubling fluctuation, with ~ C 2 sign changes 
of on-site energies at neighboring sites, is drastically suppressed. On the quantitative level, depletion of the large-i 

r / ~t 2 m ^ 2 ^<= r 2 

tail in f2(E 7 t), observed numerically in Ref. [25] can be estimated as exp ( — y^h) ~ 

[f 2 (E,t)r°, so that for 

fi c » 1 the effect is indeed dramatic. The meaning of the power (£R C ) 2 is that maintaining the period-doubling 
order requires for each site to "pay the price" to all its ~ R 2 neighbors, (iv) In 3d, the corresponding period- 
doubing fluctuation has the form V(n x ,n y ,n z ) = V(n x ) + V(n y ) + V(n z ), with V(n) — Vosign(n) exp(i7rn), so that, 
|\I>(n)| 2 oc exp [— 273(|n a | + \n y \ + \n z \)]. In lexicographic presentation [22,23] this decay manifests itself as a system 
of prominent quasi-periodic peaks with a period close to 2L 2 . In simulations [22,23] the side of the cube was small, 
L 3 = 12. Then the ALS extends over the entire system. In Fig. 2(b) we present the lexicographic plot of the analytical 
solution |^'(n)| 2 for L 3 = 12 and 73 = 0.8. We find that the shape of | \I/ (n) | 2 is remarkably close to that in numerics 
of Ref. [22,23]. 

The main message of the present paper is that, in order to test numerically, within the Anderson model, the 
predictions concerning the ALS of "continuous" theories, simulations must be carried out not too far from the band 
edges (E = ±4 in 2d and E = ±6 in 3d) where the continuous description applies. Simulations performed for E close 
to the band center reveal lattice- specific ALS that do not exist in continuous models. 

We acknowledge the support of the National Science Foundation under Grant No. DMR-0202790 and of the 
Petroleum Research Fund under Grant No. 37890-AC6. 
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FIG. 1. Shown is one quadrant of square with a side C/\/2. On-site energies in the presence of a period-doubling fluctuation 
are shown in the units of Vb; 7r-phase slip is shown only in the x-direction. Inset: solid line is a surface E(k) — E. Dotted lines 
are regions in fc-space perturbed by the fluctuation. 




FIG. 2. (a) Numerical coefficient Ce in the dependence |ln/2(-E, i)| = CEln 2 (i/72) is plotted from Eq. (11) versus the 
disorder strength, W, at different energies. Dotted line is a weak-disorder asymptotics, Ce oc W~ 2 . Inset: Ce{W) at E = 2 
is shown in the domain of the weak disorder, (b) Wave function of an ALS with a decrement 73 = 0.7 in a cube with a side 
L — 12 is shown in lexicographic order m — » L 2 (m x — 1) + L(m y — 1) + m z . 
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